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We present photonic transport in D=3 dimensional granular dissipative semiconductor media. 

Our self-consistent diagrammatic approach of the Bethe-Salpeter equation including the Cooperon 
predicts a spatially dependent correlation volume within the disordered medium that is equivalent to 
the laser’s mode extent inside the sample. Our results are relevant for the construction of disordered 
complex semiconductor micro-lasers and active Mie sphere samples. They define concrete critical 
parameters for the phase transition to lasing in disordered non-linear random media. 
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Light transport in disordered nonlinear random me¬ 
dia is a fascinating subject that culminates in research of 
random lasers [H-Q , DFB lasers @ and nonlinear optical 
fibres [ 3 ] . If these systems can be operated spectrally and 
spatially controlled they are featuring a future as large 
area coherent light sources ahead of all state of the art 
LEDs. Theoretically D=3 dimensional random lasers in 
finite dimensions confront scientists usually with tremen¬ 
dous computational efforts. Within extended granular 
systems this can only be reached by novel approaches of 
quantum field theory Q. Expanding them to photonics 
[9| and to light transport in matter under strong mag¬ 
netic fields |10l | is the aim of high-performance numerical 
physics. In this letter we develop the first theory that is 
able to describe the modal behavior of random lasers. 
The electrodynamics for transport in non-linear media is 
described basically by the wave equation for Kerr-Media 
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where the polarizability P{E) may be decomposed in 
linear and non-linear part P{E) = eo(e — 1)E -|- Pnl- 
The electric displacement inside a medium is written as 
D{E) = eoeE-\- Pnl- Kerr media are or higher order 
processes in their dependency to the electrical field E. 

P = eoX^^^E + eox^^'>E-E+... ( 2 ) 

We describe the formalism yielding the extended Bethe- 
Salpeter equation Fig[T] and lasing in steps (I) the in¬ 
teraction of the light wave with the single independent 
scatterer, (II) the propagation of the light wave in disor¬ 
dered random media, (III) the propagation with respect 
to the intensity correlation and interference effects and 
(IV) lasing. 

The propagating wave is described by the single-particle 
Green’s function Eq.Q solving the non-linear wave equa¬ 
tion. 
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FIG. 1: (a) Schematic representation of the Bethe-Salpeter 
equation, (b) Irreducible vertex including all orders of most 
crossed diagrams (Cooperon) which represents all quantum- 
coherent interference contributions in presence of disorder. 


The granularity of the medium (I) as well as it’s random¬ 
ness are implemented as spatially dependent potential of 
the permittivity €s(r). The order of non-linearity is a 
matter of the numerical selfconsistency as the indepen¬ 
dent semiconductor scatterers ll| are well described by 
the T-matrix l^, 131 for spherical Mie particles including 
the self-consistent complex refractive index. 

We use in presence a model of independent Mie scatter¬ 
ers Figl2]for the self-energy contribution S‘2. Originally 
TB| the Mie scattering coefficients of n-th order read 


^ ^ '4>n{my)A>'^{y) - mA/rviyW^^jmy) 

A/n{my)^'^{y) - m^n{y)'i>n{my)- 
In this notation m = — is the relative refractive index 

rih 

between scatterer and background, y = is the size 

parameter depending on the scatterers’ radius rs as well 
as the transport wavelength A. A>n, are Riccati-Bessel 
functions. The refractive index ns(‘I*) inherits a higher- 
order non-linearity as consequence of the selfconsistency 
Imes($). The independent scatterer approach works very 
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FIG. 2: (Left) Sketch of a disordered semi-conductor random 
laser slab open in 2 —direction. Photonic transport processes 
(red) and their time reversal (green) and internal resonances 
within independent Mie scatterers. (Right) 4-level laser rate 
equation scheme. Sinuous lines are electronic procedures, 2- 
photon pumping (green), spontaneous 732 and spontaneous 
emission 7sp, stimulated emission 721, hui represents lasing. 


well in dense random media of filling of 35% — 55%. The 
inter-scatterer dynamics is effectively mapped on the dy¬ 
namical ns($) and non-linear response even though we 
use the far field description of Mie scattering. 

The wave field solves (II) the wave equation Eq. 
(HD and builds up the intensity field {Q, 11) in 
the sample. The four point correlator $^-,((5,11) = 
r, (^)Gr r, (w)) expressed in terms of retarded 

' k+ k- k_^ '' 

and advanced Green’s function denoting the field 

and its time-reversal. We utilize the transformation 
of coordinates in center-of-motion (Q, 17) and relative 
(fc, w) momenta and frequencies 1^ k± = k ±Q/2 and 
u!± = w ± 17/2. Considering the slab geometry extended 
within the {x, ?/)-plane and finite in z-direction, the full 
Fourier transform as indicated in infinite samples 17| is 
replaced by a partial Fourier transform following the ar¬ 
gument of scale separation for field 'h and intensity $. 
'h is characterized by the wavelength A whereas change 
of the light intensity $ is characterized with the trans¬ 
port mean free path ^ derived in due course of the paper. 
In (x, 2 /)-plane the standard Fourier transform is used, 
in the limited z-direction we Fourier transform the rel¬ 
ative coordinate but the center-of-motion coordinate Z 
remains as is. 


Applying the same arguments (III) to the equation of 
motion for the intensity correlation, known as Bethe- 
Salpeter equation = G^G^[1 + f we 

obtain the Boltzmann or kinetic equation for trans¬ 
port Eq.([ 6 ]), where we introduced the following abbre¬ 
viations, also for later use, AS = E^(w_) - E«(w+), 
□E = E"^(w_) -I- S'^(w+), and equivalent expressions for 
and The term 7 ^Z", Qn, 17) is the 

irreducible vertex representing coherent photonic inter¬ 
action in the disordered granular non-linear system. The 
Ward identity, here for laser active, complex refractive 
Mie scatterers, ensures local energy conservation [l^ . 
The solution to Eq. dH]) , is derived to be the energy density 
response in its diffusion pole structure Eq.dS]) 
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FIG. 3: Coherence volume within a D=3 dimensional random 
laser slab of ZnO particles. Scatterers’ radius is r = 212 nm, 
the real part of permittivity Recs = 4.0164, the transport 
wavelength is A = 385 nm. The strongly scattering medium 
of monodisperse ZnO Mie spheres is open and lossy in the 
z-direction. Pump rate P = 1 (red), 3 (blue), 5 (yellow) in 
units of 1 / 721 . The stationary state marks the laser threshold. 
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FIG. 4: Left panel: Indeptli dependent diffusion constant 
D{Q, = 0;^). Right panel: Indepth dependent scattering 
mean free path R = - , ^ Both quantities show 

differing behavior at the samples borders ±0.5d. 


We define the correlation length ^ depending on the 
energy density $££ itself 
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Here the full diffusion constant = 0; Z) is de¬ 

rived according to D[Vl) = Dq°* — D{yL)M{VL) where 
bare diffusion Dq = f • Q)^AG is comple¬ 

mented by terms originating from the active medium 
as the scatterers Dg = ^r^AgT^Do or the background 
Db = j{uJT)^AebDo. So the diffusion constant without 
memory effects reads Dq°* = Dq + Dg + Dt- The last 
term —T^D{rL)M{Vi) contains the memory kernel M(0), 
Eq.©, that bares the Cooperon contribution and con¬ 
sequentially all interferences. Dq equals Dq where the 
imaginary part AG replaced by the real part DG. The 
renormalized density of states LDOS is derived to be 
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We use the notation and the abbreviations of [l6l| . 
A{uj) = BgAg + irgdaAgiLl) with Ag = 2[ueReGo -I- 

dfi is the differen- 
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tial resulting from the expansion in Q and fl. This 
includes the momentum integrated AG to be AGq = 
/ ^AG-(g, Z) and Ae = e, - Ug = 

and Vg = ImAe/ReAe and ^ 


ness = ^Imcf,. ve = i+A(..) 

transport velocity, whereas cp = Re 


Im(AeG^) 
Reef,, for complete- 
^ is the energy 
equals the 


phase velocity, each selfconsistently derived. The time 
scale T is discussed in detail in |16|. 

= -Wxf - Cl [d^<^>gg{Q, H)) + C2. (10) 



FIG. 5: Selfconsistent coherence length ^ for various thick¬ 
nesses d of the random laser slab. Towards strong pumping P 
a transition to a power law behavior of ^ with is found. 

Results are for identical Cs = 212 nm of ZnO (n = 2.0041) 
and G (diamond, n = 2.4). Diamond at transport wavelength 
A = 385 nm is significantly closer to the Mie resonance. 


By solving (IV) the renormalized diffusion equation 
Eq. (fTUl) coupled to the 4-level laser rate equations [1^ [l^ 
the coefficients Ci and C 2 are selfconsistently derived, and 
we arrive the spatial distribution of energy density: 
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The nonlinear self-consistent microscopic random laser 
gain 721^2 incorporates the influence of the boundary 
renormalized length scale Xd, 72 i is the transition rate 
from laser level 2 to level 1, n 2 is the electronic occupation 
number of level 2. This yields the inversion condition 
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in stationary state. Xd is the length scale implicated by 
dissipation in the bulk alone. The modification of the 
boundary specifies the relation between lasing emission 
and amplified spontaneous contributions (ASE), Eq. (ITTI) , 
relevant for light emitting diodes (LEDs). 

As results, we find a large set of selfconsistent scales, 
e.g. the diffusion constant D the scattering mean free 
path Ig the transport mean free path or coherence length 
the energy transport velocity up as well as the phase 
velocity Cp, the energy and intensity dependent LDOS 
the intensity correlator including interference influ¬ 
ences $££ as well as the laser thresholds and the gain in 
stationary state. Our results for the energy transport ve¬ 
locity Vg and the phase velocity Cp in connection to the 
LDOS Eq.(l9]) are discussed in following up work with 
respect to sub- and hyper-diffusion. 

As the central result in this paper we derive the D=3 di¬ 
mensional pump intensity dependent coherence volume 
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Fig. [3] resulting from the coherence scale ^ which is dis¬ 
played dependent to the samples depth. The inherent dif¬ 
ferential in Eg. (1101) is technically extremely pronounced 
in the decrease of coherence in Fig. [31 This result is novel 
in theory. First it is obvious that the mode loses coher¬ 
ence near the samples open surfaces (z-direction) but it is 
striking that a huge coherent transport mean free path ^ 
can be built up in the samples center of about 16, 96 /rm 
which is about factor 5 of the scattering mean free path 
and about 45-times larger than the wavelength of 385 nm. 
This length is a direct measure of the Cooperon contri¬ 
bution, which is the fully coherent photon and its time 
reversal procedure in the intensity correlator $££ that 
acts as the stimulation process of coherent emission. This 
contribution suffers due to boundary effects that are po¬ 
sition dependent whereas the monodisperse scatterers’ 
nonlinearity decreases its amount homogeneously. It is 
emphasized that coherent emission as in the laser in sta¬ 
tionary state is different from purely multiple scattered 
light, even though including nonlinearities [2l|b Conse¬ 
quentially the transport mean free path ^ in |17| is for¬ 
mally equivalent to the here derived correlation length 
however its physical meaning is fundamentally differ¬ 
ent. This is actually seen in the experiment by measuring 
second order (temporal) coherence in a Michelson exper¬ 
iment (iij i. 

The second main result of this paper is, that in sta¬ 
tionary state coherence decreases with increasing pump 
strength P whereas for comparison the selfconsistent dif¬ 
fusion constant D and the scattering mean free path Is, 
Fig. a express an inverse behavior. While diffusion is 
increased, localization (Cooperon contribution) is defi¬ 
nitely decreasing, and the scattering mean free path Is is 
large for enhanced pumping P, near the samples bound¬ 
aries [ 3 , [13 but even more in the bulk of the sample. We 
expect these findings to be measured in close up experi¬ 
ments soon. Good candidates are quasi D = 1 or D = 2 
dimensional disordered photonic crystals 3-261 disorder 
is fairly controlled and measurement access is given from 
distinct directions correlated to pure transport or out- 
of-plane scattering and emission. Also spatially varying 
gain is realizable in a controlled way . 

Qualitatively equal but inverted is the effect definitely 
present on a scale much smaller than the wavelength A 
in the depth-dependent behavior of the diffusion constant 
D as well as the scattering mean free path Ig Fig. |4l This 
conhrms hrst estimations of boundaries in earlier work by 
Akkermans et al. [3 on Anderson localization. However 
the result in Fig. |3| can not be derived if the difference 
due to boundaries would be neglected. We emphasize 
that homogeneous loss as it has been considered in 
will certainly not yield comparable results. 

In Fig. [5] we discuss the coherence length ^ and its 
transition from linear behavior to power law for identi¬ 
cal particles of the refractive index of ZnO and C (dia¬ 
mond) with scatterer radius = 212 nm. The particle 


size is comparable to the wavelenth A = 385 nm. The 
results for diamond are closer to the Mie resonant case. 
Consequentially the same pumping results in higher laser 
efficiency and an increasing nph, but also in a larger co¬ 
herence length ^ which is found in Fig. |S] The impor¬ 
tance of the Mie resonator [3 in as such is underlined 
by these findings, even though we are still far off the ac¬ 
tual resonant case. In recent work on non-equilibrium 
dynamics in ZnO Mie resonators the extreme influence 
of the resonance on the life times of electronic states, op¬ 
tical conductivity and consequentially gain is shown |32j| . 
In our considerations here, we estimate quasi-stationary 
pumping, so it is expected that for the ultrafast regime 
the effect will play even a bigger role for the dwell times 
[33 and the laser thresholds. 

Our results of numerical fit-parameter free Vollhardt- 
Wolfle theory suggest that boundaries significantly renor¬ 
malize the modes of random lasers, specifically indepth 
of the random laser sample. A mode volume is created 
on average which will result in an equivalent spectrum in 
the measurement. This volume is a characteristic length 
scale for the radom laser construction. Furthermore we 
find the inverse behavior between the characteristic quan¬ 
tities D, Is on the one hand and the coherence length ^ on 
the other. We built our framework in this paper on the 
scalar non-linear wave-equation, which is a fairly good 
approximation for particles going towards the Rayleigh 
case in homogeneous randomized samples. It will be 
subject of future work to discuss the influence of gain 
and absorption with respect to sub- and hyper-diffusion. 
With respect to quasi-order and meta glasses of complex 
resonant particles and coated spheres it will be inter¬ 
esting to investigate the selfconsistent relation between 
non-linearity and the optical angular momentum (0AM) 
for light transport as well as random lasing in station¬ 
ary state. The full vector wave description leads to ma¬ 
trix Green’s functions and requires a more sophisticated 
Ward identity to ensure energy conservation. 
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